Self-sustained current oscillations in spin-blockaded quantum dots 
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Self- sustained current oscillation observed in spin-blockaded double quantum dots is explained 
as a consequence of periodic motion of dynamically polarized nuclear spins (along a limit cycle) 
in the spin-blockaded regime under an external magnetic field and a spin-transfer torque. It is 
shown, based on the Landau-Lifshtz-Gilbert equation, that a sequence of semistable limit cycle, 
Hopf and homoclinic bifurcations occurs as the external field is tuned. The divergent period near 
the homoclinic bifurcation explains well why the period in the experiment is so long and varies by 
many orders of magnitudes. 

PACS numbers: 73.63.Kv, 72.25.Rb, 76.60.-k, 05.45.-a 



Quantum dots, also known as artificial atoms, have 
many properties of natural atoms such as discrete energy 
levels and shell structures [2, 0] ■ The physics involved in 
quantum dots is very rich because of tunability and com- 
parability of three energy scales: level spacing, Coulomb 
interaction, and thermal energy. Unlike natural atoms, 
quantum dots allow transport measurements. Many in- 
teresting transport phenomena like resonant tunnelling, 
Coulomb blockade, spin-blockade, Kondo effects, quan- 
tum conductance etc. have been observed and explained. 
Applications in nano-electronics, spintronics and quan- 
tum computing due to possible long coherence time have 
been proposed and in some cases implemented. It is 
known that both the nuclear and electron spin degrees 
of freedom in semiconductor nanostructures can be ma- 
nipulated by the hyperfine interactions In the en- 
deavor of using hyperfine interaction to manipulate the 
electron transport in quantum dots, one long-term unex- 
plained intriguing phenomenon was observed by Ono and 
Tarucha in 2004 [6| in spin-blockaded double quantum 
dots: tunneling current I under a DC-bias, as artistically 
shown in Fig. [1] oscillates with a period up to minutes 
under certain conditions. In this letter, we show that 
dynamically polarized nuclear spins (DPNS), governed 
by the generalized Landau-Lifshitz-Gilbert equation, can 
oscillate with time under a DC-bias in a magnetic field 
window. The back effect of the DPNS oscillation on the 
electron tunneling leads to the experimentally observed 
self-sustained current oscillation. 

A coherent theory should be consistent with the follow- 
ing experimental findings: 1) The oscillations observed 
only in the spin-blockaded regime in a magnetic field 
window accompany by a current jump. 2) Both the pe- 
riod and amplitude increase with the field before the os- 
cillation disappears. 3) The oscillation period can be 
tuned by the magnetic field by several orders of mag- 
nitude from seconds (the instrumental limit in the ex- 
periment) up to several minutes. 4) The oscillation is 
closely related to the motion of dynamically polarized 
nuclear spins (DPNS) [iS]. It is clear that the oscil- 




FIG. 1: (Color online) (a) Schematic setup of spin-blockaded 
double quantum dots device, (b) Four possible configurations 
in spin-blockaded regime. One spin-up electron is in the right 
quantum dot below the chemical potentials in configuration 
A. In configuration B, one electron on each of two dots form a 
spin singlet state. Configuration C is the spin singlet state of 
two electrons in the right dot. Spin-blockaded configuration 
D is a spin triplet state of two dots with one electron on 
each, (c) Schematic of field-dependence of triplet and singlet 
states. Triplet state |T _1 > cross with singlet state \S > near 
the resonance field B r (red circle). 



lation theory for semiconductor superlattices [ll[ is not 
applicable because the period would be order of 100ns 
electron tunneling time, instead of observed seconds and 
minutes. Thermal and impurity effects were also ruled 
out [6|. Although the nuclear spins in a quantum dot 
is known to have a very long (up to seconds) relaxation 
time, the relaxation process is not a periodic motion so 
that it cannot be the cause of the oscillation. 

In order to construct a sensible model, let us examine 
the plausible microscopic process of electrons and nuclei 
in spin-blockaded double quantum dots. In the experi- 
ment f||, two vertical disk-like quantum dots (InGaAs- 



AlGaAs multilayer structure) are weakly connected in 
series between source and drain (illustrated in Fig. [Th^. 
Four possible configurations in the spin-blockaded regime 
are showed in Fig. [TJj: One electron is trapped in the 
right dot (configuration A). The second electron, hop- 
ping from the source lead to the left dot, forms either spin 
singlet (B) or triplet (D) states. Electron tunneling cycle 
A ^ B — > C — > A is allowed while tunneling is block- 
aded in D, here C is spin singlet state of two electrons on 
the right dot. As shown in Fig. [TJi, three triplet states 
(\T° >, T ±:L >) are degenerated in the absence of an ex- 
ternal magnetic field, and are below the spin singlet state 
\S > by an energy J [l2| , order of inter-dot exchange en- 
ergy. Due to the spin blockade and weak coupling of dots 
with two leads, a leakage current of order of lpA, corre- 
sponding to 100ns electron tunneling time, exists. An 
external field lifts the triplet degeneracy, and \S > and 
> states shall anti-cross around certain field. The 
current experiences a jump near the crossing field B r be- 
cause spin blockade is partially removed by spin flipping 
there, and transition from configuration D to B becomes 
possible. It is known that this spin-flip process dynami- 
cally polarizes the nuclear spins, resulting in DPNS with 
a magnetization M [ijj]. This spin- flip process mediates 
also an effective spin transfer from electron spins to the 
nuclear spins. Adopting the view that dynamical motion 
of DPNS is responsible to the observed self-sustained cur- 
rent oscillation, we concentrate on the DPNS dynamics 
under the influence of both magnetic field and the above 
spin-flip process mediated spin transfer torque. 

It is well-known that the Landau-Lifshitz-Gilbert equa- 
tion governs the generic dynamics of a macro-spin 
preserving its magnitude while the so-called Landau- 
Lifshitze-Bloch equation is for the dynamics of a macro- 
spin whose magnitude can also vary. For simplicity and 
our belief that magnitude change of DPNS magnetiza- 
tion is irrelevant to the observed oscillation, we assume 
following dynamics for rh = M/\M\ fl4j |. 



dm f± _ dm _ 

= — m x H e ff + am x — h am x (m x x). 



dt " CJJ ' dt ' ~ v (1) 

Here t is in the units of (7M)" 1 , order of submicron 
second for GaAs with M = 10 6 A/m and the nuclear gy- 
romagnetic ratio 7 = 10(^4 • s/m) _1 [15j . The first term 
on the right-hand side of Eq. (TTJ) describes the Larmor 
precession around the effective field H e ff = H — Dm 
from both external magnetic field H along x-axis, and 
demagnetic field Dm (in the units of |M|). Cylindrical 
disk-like dots in experiment [f| have diagonal demagneti- 
zation factors with D x = D y = and D z = 1. In reality, 
a small difference between D x and D y exist either due 
to the inhomogeneity or deviation of dots from the per- 
fect cylindrical shape, our numerical results show that 
the physics reported here remain the same for D x ^ D v . 
The second term is the phenomenological Gilbert damp- 
ing with a dimensionless constant a. The third term 




FIG. 2: (Color online) Change of attractors as an exter- 
nal magnetic field varies: Only upper sphere was showed for 
clarity. Solid (dash) black line L+ denotes stable (un- 

stable) limit-cycles (LCs). P+ (P-) labels stable (unstable) 
fixed points (FPs), and red dot SP is the saddle point at 
rh = ( — 1,0,0). (a) For B <C B r and the only stable attrac- 
tor is P+; (b) B = 0.40T, a slightly over the semistable LC 
bifurcation where L+ and L_ are generated; (c) B = 0.42T, 
slightly below subcritical Hopf bifurcation field at which L_ 
merges with P + and becomes an unstable FP; (d) B = 0.62T, 
slightly below the homoclinic bifurcation field at which L+ 
touches SP. 



is the Slonczewski torque (per spin) along x-direction 
[ill originated from the transition from > to \S > 
mentioned earlier (Fig. [TJ:). The Slonczewski coefficient 
a = rjW/N is proportional to transition rate W from 
> to \S > and inversely proportional to the polar- 
ized nuclear number N ~ 10 5 — 10 6 Dimensionless 
coefficient t] = S±/h measures the average spin angu- 
lar moment quanta transferred from one electron to nu- 
clei. W in the Fermi golden rule approximation reads 

i in 113 



W = Wq 



ae 2 + r 2 



(2) 



Where 



1. 



AE > 
AE < 



AE = \g^sB\ — J is the level spacing between |T _1 > 
and \S > states. The effective Lande g- factor for GaAs 
is g = —0.44 [19]. fcs is the Boltzmann constant and 
T = 1.8K is the experimental temperature. The level 
broadening T of state |T _1 > is order of phonon energy of 
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lieV [HI], and W ~ 10 3 - 10 4 [201 is the typical resonant 
spin flip-flop rate. 

In the absence of the external field, all points on the 
equator of m-sphere are marginal stable fixed points 
(FPs) while m = (0,0, ±1) are unstable. Under a very 
small field, m = (—1, 0, 0) becomes the only saddle point 
(SP in Fig. [5]) and m = (1, 0, 0) the only stable attractor 
(P + in Fig. [2]). The two unstable FPs (upper P_ in Fig. 
[5]) move toward SP along the big cycle in x — z plane as 
B increases (indicated by the arrow in Fig. In terms 
of energy, a-term is always an energy sinker while a-term 
could be both energy sinker and source [3] ■ In the vicin- 
ity of P + , a-term serves as an energy source that can be 
seen from the fact that a-term tends to push the system 
away from P + . As B approaching B r , a-term may be- 
come large enough to destabilize P+. This is confirmed 
by the standard stability analysis [2lj by calculating the 
Liapunov exponent at P+. The Liapunov exponents at 
P+ become positive at about B = 0.42T, and the system 
has no stable FPs at this point. This entails the exis- 
tence of limit-cycle(s) (LCs) in two dimensions (current 
case). Indeed, our numerical calculations support this 
scenario. Fig. [5] shows the locations of various attractors 
of Eq. [I] at various B. For a field much smaller than B r 
(Fig. [2ji), all phase flows end at P+, the only stable FP 
of the system. At a critical field Bo, slightly smaller than 
0.40T, the system undergoes a semistable LC bifurcation 
(shown in Fig. (2b at which a pair of LCs, one stable (solid 
black line) and the other unstable (dashed line), appear 
simultaneously. The two LCs move in opposite direction 
as the field increases further. At a critical value B = Bi 
around 0.42T, the unstable LC merges with P + , under- 
going a subcritical Hopf bifurcation and turning P + into 
an unstable FP P_. DPNS changes from a static state 
to a LC state, the onset of self-sustained oscillation (Fig. 

The LC touches the saddle point at another critical 
value B — B2 of about 0.62T, undergoing a homoclinic 
bifurcation (Fig. [5Ji). The corresponding period diverges 
when B goes to B2 because the phase velocity vanishes at 
the saddle point. These are exactly what were observed 
in the experiment [f|. 

In order to obtain the field-dependence of period, one 
needs to locate first the LC for a given field. This can 
be done by using the Melnicov theory 2jJ, [22[ valid for 
a nearly conserved system of small a and a. It is well 
known that the energy dissipation rate [23| 



(—a— am x x) ■ dm (3) 

v dt ' w 



L(U) 



is a function of equal energy contour L{U), thus also a 
function of energy U = \{D x m^+D y m^+D z jn 2 z ) — H-m. 
f is the dissipated energy after the system moves along 
L(U) once. Interestingly, the Melnicov function of stable 
and unstable FPs and LCs is zeros by definition. Accord- 
ing to the Melnicov theory, the LC is approximated by 



one particular L(U*) that satisfies f(U*) = 0. Fig. is 
the U-dependence of f(U) for B = 0.40T, 0.42T, 0.62T 
(from top to bottom) in the vicinity of Bq, B\ and B2 
respectively. By definition, the energy contours for en- 
ergy extremes (P± minima and maxima) are points, thus 
the Melnicov function is zero there. These correspond to 
the far left and far right nodes in the figure. The origin 
(crossing of x— and y— axis) was chosen to be the energy 
of the saddle point, /-curve have positive (negative) slops 
are stable (unstable). To see this, consider a node U* 
with §jj \u* < and a slightly deviation of U from U*, 
say U > U* , U should decrease after the system moving 
along the L(U) once because f(U) < f(U*) = 0. Thus 
the system tends to push U back to U*. For B — 0.40T, 
zero slop of /([/) at the node around U = 0.1 corresponds 
to generation of a pair of LCs, leading to a semistable 
LC bifurcation (/(£/*) = 0, f£ \u* = 0). A slight in- 
crease of B lowers f(U) curve near U* = 0.1 and the 
node splits into two. The left (right) one with negative 
(positive) slope corresponds to stable (unstable) LC. The 
stable LC L+ moves towards the SP while unstable LC 
L_ moves towards P+ as B increases. The L_ merges 
with P + around 0.42T, and turns P+ into an unstable 
FP P_ (negative slop). This is a subcritical Hopf bi- 
furcation. Further increase of B to 0.62P, L + touches 
saddle point SP and become a homoclinic loop, resulting 
in a homoclinic bifurcation. After locating the LC for a 
given field in the window of [B\ , B 2 ] , the period of DPNS 
can be evaluated by 



dm 

dm 5 
ill 



(4) 



where dm/dt is given by Eq.([T]). As shown in Fig. |3b, the 
period increases monotonically with field and diverges at 
B 2 . 

The electron in triplet state (|P _1 >) plays dual roles. 
It not only glues together the nuclear spins through the 
hyperfine coupling so that nuclear spins become polar- 
ized, but also generates a torque on the DPNS. The 
self-sustained current oscillation comes from the periodic 
motion of DPNS along an LC. The LC originates from 
the instability of a static DPNS state under two compet- 
ing forces: One is energy input from the tunneling elec- 
trons that drives the DPNS away from its static state, 
a FP. The other is the dissipation of a-term that tends 
to push the DPNS to its FP. The periodical motion of 
DPNS leads to the current oscillation. This theory ex- 
plains why the current oscillation was only observed in 
the spin-blockaded regime and under an magnetic field. 
It also explains why the period can be several order of 
magnitude longer than the typical spin precession period. 
In principle, the current oscillation reported here should 
exist in systems with more than two dots in the spin- 
blockaded regime. However, it may be easier to observe 
it in a spin-blockaded double quantum dots because it is 
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FIG. 3: (Color online) (a) The Melnicov curves /({/) for 
B = 0.40T, 0.42T, and 0.62T (from top to bottom). Ze- 
ros with positive (negative) slope are stable (unstable) at- 
tractors. The far left node corresponds to (1,0,0), while the 
far right node (0, 0, ±1) and the origin (crossing of x— and 
y— axis) was chosen to be the energy of SP. (b) Field depen- 
dence of SSCO frequency cj. The period diverges at Bi- (c) 
Phase diagram in B — a/ a plane. Red, orange, and blue 
curves are for homoclinic, subcritical Hopf, and semistable 
LC bifurcations, respectively. FP, LC, LC/FP, and 2LCs de- 
note fixed-point, one limit-cycle, co-existence of limit-cycle 
and fixed-point, and two limit-cycles phases, respectively. 

easier to glue fewer nuclear spins in a smaller space. This 
is in a sharp contrast to the current oscillation observed 
in superlattices that originates from the negative differ- 
ential resistance [3], where it is only observed in super- 
lattices with more than 20 wells. It shall be interesting 
to see whether there are also field-induced Hopf and ho- 
moclinic bifurcations in Landau-Lifshitz-Bloch dynamics 
that governs spin dynamics whose magnitude can also 
vary. 

Interestingly, if a in Eq. [T] can vary independently 
from the field B, one can obtain the homoclinic bifur- 
cation, subcritical Hopf bifurcation, and semistable LC 
bifurcation curves as a function of B and a/ a in a sim- 
ilar way as what we explained earlier. Our results are 
showed in Fig. 3c. The phase diagram shows that Eq. 
Q] supports various stable phases (stable attractors), in- 
cluding stable FPs only, coexistence of a stable FP and 
a stable LC, and existence of two stable LCs. In sum- 
mary, there are a number of predictions in our theory to 
be confirmed. According to our analysis, a LC solution 
can only exist when spin transfer torque is large enough 
and self-sustained current oscillation appears and disap- 
pears when the tunneling current varies. The current 
can be controlled by the electron coupling between elec- 
trodes and dots through gates. Thus, this provides an 
experimental way of testing our theory. Also, our theory 
predicts multiple stable attractors (Fig. 3c) in certain 
parameter regions, either the coexistence of LC and FP 
or the existence of two LCs. Thus one should expect 
hysteresis loops or two oscillation periods if proper con- 
ditions are satisfied. 

In conclusion, we showed that an external magnetic 



field can induce a Hopf bifurcation at a low field and 
a homoclinc bifurcation at a high field for DPNS in a 
spin-blockaded double quantum dots, between which is 
the magnetic field window for the self-sustained current 
oscillation. The amplitude and period of the oscillation 
becomes bigger and bigger as the field increases, in good 
agreement with the experimental findings. The period di- 
verges at the homoclinc bifurcation in our model which 
explains well why the period is of several orders of mag- 
nitude larger than the fundamental time scales. 
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